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We present results from a numerical study of jV = 1 supersymmetric Yang-Mills theory using 
domain wall fermions. A set of dynamical simulations were performed for the gauge group SU (2) 
using the Wilson gauge action on 8 3 x 8 and 16 3 x 32 lattices. We considered a range of gluino 
masses (i.e., fifth dimension extents L s = 16 — 28 and input gluino mass values m/ = 0.01 — 
0.04) in order to perform chiral limit extrapolations of physical quantities. In these proceedings, 
we summarize our findings from a study of the Dirac spectrum and present new results for the 
topological charge on j3 = 2.3, 2.353 and 2.4 ensembles. 
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1. Introduction 

Supersymmetry (SUSY) is an important component of numerous proposals for beyond the 
standard model physics. Theories that possess SUSY exhibit a variety of fascinating properties 
which are of theoretical interest in their own right. Over the past two decades, a substantial effort 
has been devoted to the theoretical challenge of constructing nonperturbative descriptions of SUSY 
theories, and this has lead to enormous progress in the field of lattice SUSY. More recently, as a 
result of algorithmic and hardware advances, there has been a concerted effort devoted toward 
performing numerical simulations of such theories. Among the SUSY theories of interest-and one 
that is the focus of these proceedings-is JV = 1 supersymmetric Yang-Mills (SYM) theory. This 
theory is the only four-dimensional supersymmetric gauge theory which can be simulated on the 
lattice without numerical fine-tuning of operators. 

jV = 1 SYM consists of a gauge boson and a Majorana fermion (i.e., gluino), each of which 
transforms as an adjoint under the gauge group. The theory is confining, and possesses a variety 
of interesting phenomena which are accessible to nonperturbative study via lattice simulations. Of 
particular interest is whether or not a gluino condensate forms, which would signal the breakdown 
of a discrete chiral symmetry from Z2N — ► Z2. The presence of discrete chiral symmetry breaking 
implies the formation of domain walls, for which a known, and nontrivial relationship between the 
domain wall tension and condensate may be tested. At present, very little is also known about the 
low-lying spectrum of the theory, which is believed to consist of glue-glue, glue-gluino and gluino- 
gluino composite states [[IJ, Finally, lattice simulations can provide quantitative insights into 
how a soft SUSY breaking gluino mass affects such quantities as the string tension and spectrum. 

With conventional lattice discretizations of JV = 1 SYM, the only relevant SUSY violating 
radiative correction which is allowed by gauge invariance and the hypercubic symmetry of the lat- 
tice is a gluino mass term. However, because the continuum theory possesses a discrete Z2N chiral 
symmetry (which is an unbroken subgroup of the anomalous U(l) axial symmetry), the gluino 
mass term will be protected from radiative corrections, provided the chiral symmetry is realized on 
the lattice. For this reason, domain wall fermions (DWFs) are an ideal fermion discretization-and 
have the added advantage that, for JV = 1 SYM, they are free of sign problems [Q]. With the use of 
DWFs, supersymmetry emerges in the continuum and chiral limits in a controlled and theoretically 
understood fashion without fine-tuning. 

We have recently performed dynamical simulations of jV = 1 SYM using DWFs Our 
study may be, to an extent, viewed as a continuation of the early exploratory simulations of [0]. In 
our work, we performed measurements of basic quantities and established important benchmarks 
which are crucial for future studies of the theory. These include: 1) a study of the residual chiral 
symmetry breaking (i.e., the residual mass m res ) and the systematic errors associated with the chiral 
extrapolation of the gluino condensate, 2) establishing the lattice scale and 3) performing a detailed 
investigation of the Dirac spectrum. 

Here, we present some of our results from a study of the Dirac spectrum, as well as new results 
for the topological charge. The primary motivation for studying the eigenvalues of the Hermitian 
DWF Dirac operator is several fold. First, we would like to verify that the DWF formulation is 
working properly. This entails verifying that the low energy modes are indeed bound to the fifth 

A similar study of ,yV = 1 SYM using DWFs was recently reported in M. 
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Table 1: Simulation parameters, gluino condensate (qq), average plaquette (P), residual mass m res and 
Sommer scale ro for a subset of ensembles found in [Q, 0]. Roman numerals serve as ensemble identifiers. 
The Sommer scales quoted for Ensembles II-IV are estimates obtained from L s = 16, 20 and 24 ensembles. 





V x T 


/3 


L s 


m f 




(P) 






I 


8 3 


x8 


2.3 


24 


0.02 


0.006806(34) 


0.73202(32) 






II 


16 3 


x32 


2.3 


28 


0.02 


0.0063346(33) 


0.731688(32) 


0.14834(13) 


3.3(1) 


III 


16 3 


x32 


2.353 


28 


0.02 


0.0057106(61) 


0.743853(38) 


0.10269(18) 


4.3(1) 


IV 


16 3 


x32 


2.4 


28 


0.02 


0.0049179(82) 


0.752951(27) 


0.06513(17) 


5.3(1) 



dimension boundaries, and that those modes properly describe the desired chiral physics, reflected 
by the nature of the matrix elements of the physical chirality operator. Second, the eigenvalues 
provide a qualitative measure of the proximity to the continuum limit; in this limit the positive 
and negative eigenvalues of 75 p should be paired, however, this pairing is broken by the lattice 
spacing. Finally, the eigenvalues provide an independent method for determining the residual mass 
and gluino condensate, where the latter is achieved by exploiting the Banks-Casher relation [S]. 

By measuring the topological charge, we are able to monitor topological ergodicity in our 
simulations. Furthermore, on a configuration by configuration basis we may also check whether 
or not the near-zero modes of the Hermitian DWF Dirac operator are consistent with the gluonic 
definition of the topological charge. 

2. Ensembles 

We have performed a series of numerical simulations of jV = 1 SYM for the gauge group 
SU (2). Simulations were performed using DWFs |^, [T(|] and Wilson gauge action on 8 3 x 8 and 
16 3 x 32 lattices. Our current simulation parameters include fifth dimension extents ranging from 
L s = 16 — 28, and input gluino mass values ranging from nif = 0.01 — 0.04. The majority of our 
simulations were performed at a single lattice spacing, which corresponds to /3 = 2.3. However, 
several ensembles were generated at finer lattice spacings, which correspond to j8 = 2.353 and 
j8 =2.4. The ensembles considered in these proceedings are a subset of those discussed in [|], ||] 
and are listed in Table ([[]). 

In order to relate the lattice parameters to the physical system, it is important to recall that 
the lightest states in jV = 1 SYM are the QCD analogs of glueballs, the r\' and it's corresponding 
fermionic superpartner (for which there is no QCD analog). The lightest states therefore have 
masses on the order of the inverse Sommer scale (r^ ), and the ratio rr,/L provides a measure of 
the finite volume errors. For our current ensembles, we find tq/L 0.18 — 0.30 for j8 = 2.3 — 
2.4, where ro in largely insensitive to m/ and L s . On the other hand, the supersymmetric limit is 
controlled by the finite lattice spacing errors and proximity to the chiral limit. The bare gluino 
mass is given by m g = rrif + m res , and hence the chiral regime corresponds to: m g rQ << 1. For our 
ensembles, the bare gluino mass is dominated by the residual mass contribution, i.e., m res /mf ~ 
<^(10), which is strongly dependent on /3 and L s . At L s = 28 the gluino mass ranges from m g rQ 
0.35 — 0.5 for j8 = 2.4 — 2.3, suggesting that we are at most only moderately within the SUSY 
regime for our current set of ensembles. 
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3. Dirac spectrum 

We measured the lowest 60 eigenvalues A# of the Hermitian DWF Dirac operator on 8 3 x 8 
lattices using the method of Kalkreuter-Simma. The details of the algorithm are described in [Q, [Tjj ]. 

The Hermitian DWF Dirac operator is given by Dh = R^D, where D is the DWF Dirac op- 
erator and R 5 is the fifth dimension reflection operator defined in We measured the eigenvalues 
of Dh for five different values of the valence mass m v , ranging from —0.16 to —0.12. Negative 
values for the valence mass were chosen on the order of m res in order to minimize the gluino mass 
on each background gauge field configuration. The valence mass dependence of the / th ei genvalue 
A H i was then fit to the reparameterized Taylor expansion JT2| ] : 

A^iK) = n 2 5 i [A, 2 + (m v + Smi) 2 ] + 0(m 3 v ) . (3.1) 

Fig. ([[]) shows a plot of A H (m v ) and fit results for the lowest ten eigenvalues obtained from a 
representative configuration in Ensemble I. 



In the continuum limit, the Hermitian Dirac operator 75 (p + m g ) has eigenvalues ± J A 2 + ra 2 ,, 
where rb'A are eigenvalues of p and A > 0. The zero modes of p correspond to the unpaired 
eigenvalues +m or — m of the Hermitian Dirac operator. At finite lattice spacing, Dh does not 
possess ± pairing of eigenvalues, and the apparent lack of pairing in Fig. ([!]) indicates that the 
Dirac spectrum is, at least qualitatively, not very continuum-like. 



Using the functional form of Eq. (3.1 ), the chiral condensate may be expressed as [12]: 



-fe) = T^(l A;2+ 7^ m 5m;)2 )- M 



The form of Eq. (3.2) allows us to identify A,- with the four-dimensional eigenvalues of D, and 
as a parameter that characterize the residual chiral symmetry breaking effects on an eigenvalue by 
eigenvalue basis. Fig. (||) shows a scatter plot of the extracted pairs (5m/, A,-) for the lowest 60 
eigenvalues of the entire ensemble. From this plot, it is evident that there are approximately ff(\0) 
near-zero modes for the entire ensemble, and that the majority of eigenvalues satisfy A,- << m g . 
Note that in the regime A «m g , eigenstates of 75 (fl+m) are near-eigenstates of 75. This suggests 
that the low-lying eigenmodes of Dh are bound to the walls of the fifth dimension. We have 
confirmed this behavior directly by studying the four-dimensional norms of the eigenstates of Dh 
as a function of the fifth dimension 

Finally, the fit parameter 5m, characterizes the residual chiral symmetry breaking effects of 
finite L s . We therefore expect a relationship between the distribution of 5ra,- values, p(Sm), and the 
residual mass. Specifically, in [|2|] it was observed that p(8m) was highly peaked around m res for 
quenched two flavor QCD. Such behavior is not evident from Fig. (fj), indicating that the residual 
chiral symmetry breaking is perhaps too large to be characterized solely by m res . 

4. Topological charge 

We measured the topological charge on 8 3 x 8 and 16 3 x 32 ensembles using both fermionic 
and gluonic definitions of the topological charge. The former is given by 

e^ocTrD^ 1 , (4.1) 
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up to a normalization factor. In order to compare this with the gluonic definition of the topological 

(f) 

charge, we average Q top over a smoothing window of size n smoot h- The latter definition is deter- 
mined using a two step procedure. First, the gauge fields are smoothed via APE smearing in order 
to remove ultra-violet noise: 



U^x)^¥ su{2) 



;i-c_)^(x) + ^p £ U v (x)U^x + v)Ul(x + ^) 



(4.2) 



where (x) represents a gauge link spanning sites x and x + pL, c smear = 0.45 is the smearing coef- 
ficient, and the sum over v takes both positive and negative orientations. After each smearing step, 
the gauge field is projected back onto the gauge group SU (2). After a total of n smear applications 
of Eq. ( |4.2[ ) are performed, the topological charge is then computed using 

Q$ p = ^Jd A xTrFF(x), (4.3) 

where F flv is the five-loop improved (51i) definition of the field strength tensor built from lxl, 
1x2, 1x3, 2x2 and 3x3 clover-leaf terms [O]. The use of APE smearing over the 51i cooling 



method of [ ]13[ ] was adopted by [|14|] in the context of QCD after finding that the results for the 
topological charge were comparable. For convenience, we chose to use APE smearing as well, 
however, in the future we intend to check our analysis with cooled lattices. 

Fig. (||) shows a plot of the topological charge as a function of n smear for some represen- 
tative configurations in Ensemble IV. Beyond n smear ^ 100, the topological charge takes integer 
values, and infrequently fluctuates between integers. These fluctuations are likely attributed to the 
disappearance of small instantons that are unstable under the APE smearing procedure. Fig. (|J) 
shows a histogram of the topological charge for three different values of the coupling (Ensembles 
II-IV). The Gaussian shape of the distribution indicates that the topological charge is adequately 
sampled in our numerical simulations. A comparison of the fermionic and gluonic definitions of 
the topological charge is provided in Fig. (||) for all three coupling values. The two definitions 
yield consistent results at the weakest coupling (Ensemble IV), however at stronger couplings, the 
disagreement between definitions becomes increasingly evident. 

On a given ensemble, we may directly compare the near-zero modes of D# with the gluonic 
definition of the topological charge. Fig. @ is plot of the topological charge as a function of 
trajectory number for the same ensemble as Fig. (^), allowing for a direct test of the Atiyah-Singer 
index theorem. Strictly speaking, at finite L s , the DWF Dirac operator does not obey an index 
theorem, however, for sufficiently large L s , one expects that it will do so approximately. Based on 
the relatively few near-zero modes in Fig. @ and the frequent sampling of non-zero topological 
charge sectors in Fig. @, however, there is little evidence that we are in such a regime in the case 
of our strongest coupling ensemble. 

It would be interesting to test whether or not the index theorem is realized on our weakest 
coupling ensembles, since for these the comparison of fermionic and gluonic definitions of the 
topological charge in Fig. (|6|) exhibit greater consistency. This consistency is a likely indication 
that the low-lying spectrum of D# has the desired continuum pairing of eigenvalues and that zero- 
modes are consistent with the topological charge. However, because the adjoint fermion Dirac 
operator has two-fold degeneracy due to a charge conjugation symmetry, measuring eigenvalues 
on our 16 3 x 32 lattices is too costly with use of our current eigenvalue measurement code. 
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5. Conclusion and outlook 

Our Dirac spectrum and topological charge results are all consistent with the existence of large 
residual chiral symmetry breaking attributed to finite L s . The results underscore the necessity for 
reducing the residual chiral symmetry breaking by going to weaker coupling, increasing L s , and/or 
exploring improved lattice actions. We are optimistic that with a combined approach which utilizes 
all three techniques, a successful study of jY = 1 SYM near the SUSY limit can be achieved. 
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Figure 1: Lowest 10 eigenvalues A# as a function 
of m v for a typical configuration in Ensemble I. Solid 
curves represent fits to data using Eq. (3.1). 
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Figure 2: Scatter plot of the fitted pairs (A, 8m) ob- 
tained from the lowest 60 eigenvalues of all configu- 
rations in Ensemble I. 
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Figure 3: Topological charge (gluonic) as a function 
of n S mear for 10 different configurations in Ensemble 
IV. 



Figure 4: Histogram of the topological charge (glu- 
onic) evaluated at n smear = 60 for Ensembles II, III 
and IV. 
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Figure 5: Comparison of gluonic (n smear = 60) and 
fermionic (n smoot f, — 5) definitions of the topological 
charge for Ensembles II, III and IV. 
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Figure 6: Topological charge (gluonic) evaluated at 
nsmear = 60 as a function of trajectory number for 
Ensemble I. 
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